Inside looking out: probing JIMWLK with BFKL calculations. 
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We investigate the relation between the eigenvalues and eigenfunctions of the BFKL and 
JIMWLK/KLWMIJ Hamiltonians. We show that the eigenvalues of the BFKL Hamiltonians are 
also exact eigenvalues of the KLWMIJ (and JIMWLK) Hamiltonian, albeit corresponding to possibly 
non normalizable eigenfunctions. The question whether a given eigenfunction of BFKL corresponds 
to a normalizable eigenfunction of KLWMIJ is rather complicated, except in some obvious cases, 
and requires independent investigation. As an example to illustrate this relation we concentrate on 
the color octet exchange in the framework of KLWMIJ Hamiltonian. We show that it corresponds 
to the reggeized gluon exchange of BFKL, and find first correction to the BFKL wave function, 
which has the meaning of the impact factor for shadowing correction to the reggeized gluon. We 
also show that the bootstrap condition in the KLWMIJ framework is satisfied automatically and 
does not carry any additional information to that contained in the second quantized structure of the 
KLWMIJ Hamiltonian. This is an example of how the bootstrap condition inherent in the t-channel 
unitarity, arises in the s-channel picture. 



INTRODUCTION 

In this paper we study the relation between the BFKL perturbative resummation of leading logarithms [l[ at high 
energy and the JIMWLK/KLWMIJ evolution 0,0] equations, which take into account the physics of saturation 
and multiple scatterings [4] . With a slight abuse of language we call BFKL the second quantized formulation which 
leads not only to the BFKL equation for two gluon exchange, but also to the whole set of BKP equations for t- 
channel exchanges with an arbitrary number of gluons [5] . The basic correspondence between the elements of the two 
approaches was established in @ and we will review it briefly below. 

Here we are interested in the question how to relate the eigenvalues and the eigenfunctions of the BFKL Hamil- 
tonian and their counterparts in the JIMWLK/KLWMIJ theory. We will show that the eigenfunctions of the 
JIMWLK/KLWMIJ Hamiltonians when expanded in Taylor series in the appropriate variable (p for JIMWLK and 
S/5p for KLWMIJ) are eigenfunctions of the BFKL Hamiltonian. However most of the BFKL eigenfunctions when 
resummed into solutions of JIMWLK/KLWMIJ are not normalizable. The question which are, and which are not, 
cannot be settled within the BFKL framework. We also discuss how to calculate higher corrections to the BFKL 
eigenfunctions. As an example we concentrate on the reggeized gluon Q. We calculate corrections to the reggeized 
gluon wave function. We show that terms with up to two gluons in the t-channel reggeize due to the bootstrap 
condition inherent in this approach. We also point out that in the JIMWLK/KLWMIJ framework the bootstrap is a 
necessary consequence of the hermiticity of the JIMWLK/KLWMIJ Hamiltonian, and does not require a specific form 
of the emission kernel. We calculate corrections to the eigenfunction beyond the two gluon exchange approximation 
and thereby find the screening correction which appears when at least three gluons are exchanged in the t-channel. 
This is an interesting example of the interrelation between the t-channel unitarity, which is the origin of Reggeons 
and the s-channel one which is the inherent feature of the JIMWLK/KLWMIJ approach. 

Let us start by recapitulating the JIMWLK/KLWMIJ formalism. In this approach one considers the scattering of 
a projectile hadron on a hadronic target. The projectile is described by a distribution of color charge density in the 
transverse plane pp(x), while the target is viewed as an ensemble of the color fields (Xt{x) with probability densities 
PF p [/9p] and Wrlctr] respectively. The second quantized S matrix operator is given by its eikonal expression (see 
Fig. P 

S = exp^iJ dy~ J d 2 x p a p(x,y~)a>T(x,y~)^ (1) 
and the forward S matrix element at rapidity Y is given by the functional integral 

S Y = J Da a T J dp P W£[p P ] W T [a T ] expji^ dy' J ' d 2 'x p a P (x , y~ ") ct a T {x , y~ )| (2) 
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FIG. 1: The S matrix in the JIMWKL approach with averaging specified by Eq. p|. 



The rapidity evolution of basic physical observables (including the S'-matrix) is determined by the second quantized 
evolution Hamiltonian Hrft (the following pertains to the projectile evolution - we drop the subscript on p for 
brevity) 

- ^yW\p) = H RFT [ P , j-]W\p] (3) 

The evolution Hamiltonian Hrft depends on two unitary matrices, 

S(x) = VexpS^iJ^ dy-T a a a (x,y-)^ R(x) =7>cxpj^ daT * T a \ . (4) 

Here the field a is the projectile color field analogous to ax- It is related by a nonlinear transformation to the 
projectile color charge density p 

a a (x,x-)T a = g 2 -L(x - y) {S\y,x~) P a (y,x~)T a S{y, X -)} . (5) 

with ^r(x,y) = ± \n[(x - y) 2 p 2 ] and 

S(x,x-) = Vexpli [ dy-T a a a (x,y-)\ (6) 



The S- matrix at arbitrary rapidity can therefore be represented in terms of the eigenvalues of Hrft- Its evolution 
is given by 



dSy 



J Da^ J dp P H RFT [p,-^-}W^[p}W T [a T }exp^i J dy~ J d 2 x p a (x,y-) a a T {x,y-) 



(7) 



dY 

Given the set of complete eigenfunctionals of Hrft 

Hrft^M = WiVi\p] (8) 
one can expand the probability distribution at the initial rapidity Y Q as 

wv„[p] = I>*<M ( 9 ) 

i 

and therefore 

Sy =J2e- Ui(Y - Yo HPi (10) 
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with 



ft 



J Da a T J dp^i[p] W T [a T ) ex p{*^ d V~ J ^ x P a 



(x,y )a%(x,y ) 



(11) 



Given that the Hamiltonian H R ft is positive definite [Io| , the asymptotic behavior of the S- matrix is governed by 
the lowest lying eigenvalues Wj. The study of the spectrum of H R ft is therefore of direct relevance to understanding 
the high energy behavior of physical amplitudes. 

Two limits of Hamiltonian H R ft have been widely discussed in the literature. One is valid in the limit of dense 
projectile, but allows only exchange of two gluons with the target 0. In the other limit one assumes that the projectile 
is dilute, but resumms all possible multiple interactions of the partons of the projectile with the target [3j . The two 
limiting cases are related by the dense-dilute duality transformation 11| and have therefore identical spectra. In the 
rest of this paper we choose to study the KLWMIJ Hamiltonian. 



H 



KLWMIJ 



2tt 2 



K xyz {J a L {x)Jl{y) + J a R (x)J a R (y) - 2J L (x)R a z b J b R (y)} 



x.y.z 



with the kernel 



_ (x-z)j(y-z)j 



and the left and right rotation generators 



Jftx) = -tr 
J a R {x) = -tr 



SRl 
R x T a - 



-T a R, 



5Rl 



Alternatively one can write 



Hklwmu = £5 [ K xyz \j^(x)J^(y) - 2J a L (x) (R z - l) ab J b R (y)\ 

-» Jx,y,z 1 " ' 



with 



J^{x) = J a R {x) - J a L {x) 



(12) 

(13) 

(14) 
(15) 

(16) 
(17) 



The operators Jy{x) are the generators of SUv{N c ) - the vector subgroup of SUl(N c ) ® SU R (N C ). 

As noted above, the JIMWLK Hamiltonian is obtained from cq. (IT2l by the dense-dilute duality transformation 



R(x) -> S(x) 



(18) 



We note that a generalization of H R ft which interpolates between Hklwmu and Hjimwlk has recently been 
derived (l2| . However due to its complexity we will not deal with it in the present paper. 

The Hamiltonian Hklwmu is the limit of H R ft at low color charge density. The left and right rotation operators 
J R and Jl eq. (|14|) that appear in eq. (|12p are in fact just the color charge density in the hadronic wave function and 
its conjugate respectively [9(. 

Some aspects of the structure of the spectrum of Hklwmu were discussed in [lp| . In particular we know that 
Hklwmu is positive definite since it can be written as 



Hklwmu = £5 J Q- \z)Q a t (z) 



with the Hermitian amplitude 

It has two states with zero eigenvalue: 



% [R« b (z)-R"»(x)]J R (x) 

x (x- z y 



(19) 



(20) 



H K LWMu\Yin) = 0; H K LWMu\Yang) = 



(21) 
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The state \Yang) corresponds to the physical vacuum, that is to the state which is annihilated by the color charge 
density 



J a L {x)\Yang) = J a R (x)\Yang) = 



while \Yin) represents the black disk state 



R a \x)\Ym) = 5 ab \Yin) 



(22) 
(23) 



We will find it convenient to work in the basis of eigefunctions of the operator R. Written in this basis the wave 
functionals of the two states are 



(R\Yang) = 1; (R\Yin) = 5 (R ab (x) - 5 ab ) 



(24) 



Each one of these states sustains a tower of excitations above it. The RFT states "close" to \Yang) correspond to 
physical QCD states with small number of particles in the projectile wave function. This interpretation stems from 
the fact that, as discussed in detail in Q the probability distribution W has a convenient representation 



W[ P ] = V[R]6\p] 



(25) 



and thus in the R basis W is simply a regular function of R. Every factor of R in W corresponds to a gluon in the 
projectile wave function. Thus expansion of W in powers of R — 1, is equivalent to expansion in the number of the 
projectile gluons that scatter on the target. 

In this paper we are interested in the eigenstates of Hk lwmij which have similar structure, namely in the R basis 
the eigenfunctions are regular functionals of R(x) which can be expanded in powers of R — 1. Since R is a regular 
function of S/Sp, these same states can also be expanded in powers of 5/5p. Clearly the state \Yang) is one of those 
states since the wave function in this case is simply a constant. 

On the other hand the state \Yin) does not fall into this category. Its wave function is not expandable in powers 
of R — 1. The same is true for other states "close" to it. As explained in [l(| those states correspond to "holes" in 
the black disk and their eigenfunctions are expandable in powers of p rather than S/Sp. 

The discussion of this paper pertains directly only to the \Yang) - like states. 

The derivation of the Hamiltonian Hk lwmij [2| does not assume anything about the strength of the target fields. 
If one interested in the situation when the target fields are small, one can expand Hk lwmij in powers of S/Sp. The 
expansion in powers of S/Sp is equivalent to expansion in powers of the target field, since powers of S/Sp turn into 
powers of ax in the calculation of the scattering matrix eq.([7]). The leading order of this expansion gives Hbfkl, 
which is the second quantized Hamiltonian that generates the high energy evolution in the BFKL framework. The 
form of Hbfkl is well known (see for example [13|). For completeness we present the derivation of Hbfkl in the 
appendix, carefully keeping track of the path ordering in the definition of R(x). Although this path ordering is not 
relevant in many cases [6(, in general it cannot be neglected. The result is 



Hbfkl 



2ir 2 



K xvz {T a T\ d p a x 



8 _S_ 

5p% Sp% 



S J_ 



Pi 



where 



(x,x ) 



In the "normal ordered form" this reads 
Hbfkl 

with 



a s 
~2ir 2 



[ K xyz (T a T b ) cd 

J xyz 



S 5 
Sp% 5p c z 



5 5 
Jo~l~Jp~i 



P> b y + 



: $p a z 



(3 X - Z = 



XsNc 

2vr 2 



S' 2 (x — z) I K(x, x, u) — K(x, x, z) 



a s N c 
2tt 2 



S 2 (x-z) 



{x-zf 



(26) 



(27) 



(28) 



(29) 



The question we want to address is what is the relation between the eigenvalues and eigenfunctions of Hbfkl and 
Hklwmij ■ More importantly, to what extent can we use the results of calculations performed in the framework of 
BFKL evolution to get information about the spectrum of Hklwmij- 
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FROM BFKL TO KLWMIJ 



The calculation of the spectrum of Hbfkl is equivalent to solution of the complete set of BKP equations. 
The BFKL Hamiltonian is a homogeneous function of the coordinates (S/Sp) and momenta (p) and thus its eigen- 
functions are pure powers. Taking an eigenfunction in the form 

9 A [6/Sp] = ( G£°»"°» ( Xl ... Xn )jL...-lL (30) 

Jxi...x n °Pxi °P*l 

and acting on it with Hbfkl leads to an eigenvalue equation of the form 



2tt 



\ 1 / rparpb\ 

^i^jK 1 1 Jaictj 



-XiXjZ {p^l ■ •■•^n ) H~ / K X y Xi G a (^1 • --X.-.y ...x n 

Jxy 

I -^Xiyxj G (xi . . .X{ ...y...X n ) / K$XjXi G A * • • *^n) 

J v J X 

0xB t G%"- a -" b -" an (xi...x...x n ) =u A G^- a "(x 1 ...x n ) (31) 



Eqs. (|31l) are precisely the BKP equations [5| for the rapidity evolution of n - gluon exchange in t-channel. The index 
A denotes various quantum numbers that characterize the eigenfunction Ga, in particular total momentum, color 
representation, charge conjugation, parity and so on. 

As mentioned above, expanding Hklwmij and \& is powers of S/Sp is equivalent to expanding the S matrix eq.@ 
in powers of the target color field cut- Every factor of olt represents an exchange of a gluon in i-channel between the 
projectile and the target. Thus physically, expansion in powers of 5/Sp is equivalent to expansion in the number of 
gluons exchanged in the t - channel [6|. Consequently, the quantum numbers denoted by the index A in eq. (|31|) arc 
the quantum numbers of the n - gluon state exchanged in the i-channel. 

For n = 1 the solution of eq. (|3T]) is the reggeized gluon. In this case the color representation of the exchange is 
obviously adjoint and the eigenvalues are characterized by the transverse momentum. For the singlet two gluon state, 
n = 2, this is the celebrated BFKL equation and the solution is the BFKL Pomeron. For arbitrary n in the large N c 
approximation these equations have been extensively studied in LJ] where it was shown that the spectrum of the n 
gluon state is described by an integrable spin chain. 

Although the full spectrum of eq.(|3T|) is not known, the salient features are the following. In the reggeized gluon 
sector (n — 1) the eigenvalues are nonnegative. The lowest eigenvalue is vanishing and corresponds to zero transverse 
momentum exchange t. At nonzero t the eigenvalues are logarithmically infrared divergent. For the BFKL Pomeron 
(singlet n = 2 exchange) the lowest eigenvalue is actually negative, corresponding to the growth of the amplitude 
at high energy with the famous BFKL intercept. At n > 2 the spectrum is rich. Importantly the lowest eigenvalue 
for the color singlet exchange is always negative and grows proportionally to n, corresponding to the n/2 Pomeron 
exchanges [1H 

Now let us consider Hklwmij- The eigenvalues and eigenfunctions as noted above are determined by solving 

Hklwmij* a[R]=ua*a[R] (32) 
Suppose we try to solve this by Taylor expanding ^ a in powers of S/Sp. We thus write 

* A = *°[A] + *^] + ... (33) 

To find ^ we need to act on it with Hklwmij also expanded in powers of S/Sp. It is clear from the mechanics of the 
expansion of Jl and Jr given in the appendix, that to all orders in S/Sp, both the charge densities are proportional 
to the first power of p, and only the power of S/Sp in Hklwmij grow order by order. Therefore the Hamiltonian 
Hklwmij can be written as 

Hklwmij = Hbfkl + Hi + ... (34) 
Here the first correction to the Hamiltonina, Hi schematically has the form 

^=^(£)V + A(£) 2 P (35) 
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The structure of ecis. (|33l34[) is such that in the leading order the wave function ^^[/-] must be an eigenfunction of 
Hbfkl- Thus the leading order equation determines uia completely, and the role of the higher order equations is 
only to determine the higher order Taylor series terms in the expansion of the wave function ^^[i?]. E.g. having 
found ^ A as an eigenfunction of Hbfkl with eigenvalue lua, to first order one has 

*>] = 5 H^° A (36) 

and so on. 

It thus appears that the eigenvalues of Hklwmij can be obtained exactly from the leading order approximation 
- the BFKL Hamiltonian. The catch however is that we do not know from the BFKL calculation per se whether 
the wave function corresponding to a given "eigenvalue" will turn out to be normalizable or not. Within the BFKL 
approximation itself, none of the wavefunctions are of course normalizable since they are simple monomials of 5/Sp. 
However only normalizable eigenfunctions of Hklwmij should be used in the expansion of the iS-matrix eqs. (19H0p . 

A similar situation is encountered in simple quantum mechanics. Take for example a one dimensional harmonic 
oscillator 

h = \ (P 2 + x 2 ) (37) 
Let us now try to find its ground state wavefunction in Taylor expansion. We take 

V; = 1 + ax 2 + cx 4 + ... (38) 

Acting on it with the Hamiltonian we get 

hip = -a+ ^x 2 - 6cx 2 ... (39) 



and the Schroedinger equation 



Thus for arbitrary lo we simply have 



a+ i(l - 12c)x 2 = oj(1 + ax 2 ) (40) 



a = -u; c=-^(1 + 2oj 2 ) (41) 

So there is a solution for every possible value of u> (in fact there are two solutions for each w, but the other one is 
an odd function of x and is outside our initial ansatz). This calculation however does not carry a lot of information 
since we do not know a priori which of the so found functions are normalizable when the Taylor series is summed to 
all orders. We know of course that the spectrum in fact is discreet and therefore most of the "eigenvalues" do not 
correspond to normalizable eigenfunctions. 

The situation in the KLWMIJ-BFKL system is similar in this respect. As we have noted above, Hklwmij is 
hermitian and positive definite, thus its eigenvalues have to be positive. On the other hand many of the eigenvalues 
of Hbfkl are negative (including of course, the Pomeron). We are therefore assured that those eigenvalues do not 
correspond to normalizable eigenfunctions. As for the positive eigenvalues of Hbfkl, it is tempting to surmise that 
they correspond to normalizable eigenfunctions of Hklwmij- Unfortunately, we have no right to do so. The question 
whether resummed Taylor series is normalizable or not is very complicated and can not be answered in any finite 
order of Taylor expansion. We note however that the Taylor expansion for Hklwmij is in a subtle way different 
from that for the harmonic oscillator. In the later case the leading order of the expansion puts no restrictions at all 
on possible eigenvalues. In the KLWMIJ case however, the leading order itself leads to an eigenvalue problem, so that 
not every ua is allowed. 

Even though it is not clear whether the eigenfunctions of Hbfkl give rise to normalizable eigenfunctions of 
Hklwmij , it is still interesting to illustrate the procedure discussed above by some concrete example. In the following 
section we will therefore consider higher order in S/Sp corrections to some eigenfunctions. We will concentrate on the 
eigenvalues corresponding to the reggeized gluon, which is the simplest eigenfunction of Hbfkl- 
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THE REGGEIZED GLUON AND THE BOOTSTRAP. 
The Reggeized gluon 

Let us look for the eigenstate of Hklwmij which corresponds to the quantum numbers of one gluon exchange. 
The state must belong to the adjoint representation of SUv(N c ) and its wave function when Taylor expanded should 
start with the linear term in 5/5p. Thus we take 

*°=J* x Mjjh) (42) 

Acting on it with Hbfkl we obtain the eigenvalue equation 

/ Pxz<t>z = ^4>x (43) 
J z 

This is solved by 

<P{x) = e lqx (44) 

with the eigenvalue 

"< = l;l/ k mh? < 45 » 

This is nothing but the gluon reggeization. 

To calculate first correction to the reggeized gluon state we have to consider V? 1 which is quadratic in S/Sp. We will 
perform the calculation in a slightly different way which is technically simpler. We know that the wave function ^ 
at the end of the day should depend only on R. Therefore it makes sense, rather than taking an arbitrary quadratic 
function, to choose such a function that itself can be obtained from expansion of R. Since ^° is simply represented 
as expansion of R 

%[S/Sp] = J d 2 x^(x)tr (T a i?) first ordcr (46) 

we will take a simple guess 

% [S/Sp] + ttj [S/Sp] = J d 2 x4>{x)tr (T a R) Rlst+second ordcr (47) 

and will show that it indeed satisfies the KLWMIJ eigenvalue equation to second order in S/Sp. The matrix R is 
taken here in the adjoint representation. 

Bootstrap of the antisymmetric adjoint 

We now consider the action of Hklwmij in the state 

G\ = J d 2 x4>{x)tr (T a R) (48) 

The basic elements we need is the action of the left and right rotation generators on the matrix R 

[J a L (x), R(y)] = T a R(y)S 2 (x - y); [J a R (x), R(y)] = R(y)T a S 2 (x - y) (49) 
With this it is easy to calculate the action of the real and virtual parts of H KLWMIJ : 

2Jt(x)[R{z) - lY d J d R (y)tr[T a R u ] = 2[R Z - ir d tr[T a T c R x T d ]8(x - y)5(y - u) (50) 

and 

\Jl{x) - J c R {x)][Jl{y) - J R (y)MT a R u ] = tr {T"[T C , [T c , R x ]]}5(x y)5{y u) (51) 
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Hence 



Using 



we write 



H 



KLWMIJ 



G\ = - f K uuz <f> u {2[R Z ~ l] cd tr[T d T a T c R u ]+tr{T a [T c , [T C ,R U ]}}} 



tr {T a [T c , [T c , R u }}} = N c tr[T a Ru] 



HKLWMIjG a A 



/ Kuuz<t)u { 2[Rz ~ l \ Cdtr [T d T a T c R u ] - N c tr[T a Ru]} 
= 4 / K uuzK {-2[R Z - l] cd tr[T d T a T c (R u - 1)] + N c (tr[T a R u ] - tr[T a R z ])} 

" J U,Z 

We now have to expand this to second order in S/Sp. The matrix R is expanded as 



(52) 



(53) 



(54) 



r: 



ab rab 



T, 



ab 



du 



Sp c (u,u ) 



rpaerpeb 



du 



Sp c (u,u 1 ) Sp d (u,u 2 ) 



(55) 



The crucial observation is that the first and second order terms come only from the second term in eq. (l54l) . For first 
order term this is obvious. The second order contribution from the first term in eq. (|54j) is proportional to 



T^ d tr[T d T a T c T 



Sp b u Sp e z 



However, using the properties of the adjoint SU(N) generators 



tr(T a T b ) — f aa P f ba P = N c 6 ab tr(T a T b T c ) — —i f aa/3 fP b ~< p ca — % c j abc 



one can easily show that 



T^ d tr[T d T a T c T b ] = 



and so this contribution vanishes. Thus to second order the eigenvalue equation is 



aNc 
2tt 2 



K, 



uuz^ru 



_ -)+-r 

6 pi 5ft ' - 1 



x~>y- 



s 



Sp b (u,x ) Sp c (u,y ) Sp b (z,x ) Sp c (z,y ) 



_ + _ fabc 



x~>y 



Sp b (u, x ) Sp c (u, y 



(56) 



(57) 



(58) 



(59) 



Obviously, this is satisfied as before by cj>(x) of ea. (l44l) with the eigenvalue eq. (|4"5|) . According to our earlier discussion, 
the eigenvalue determined in the leading order does not change order by order. An interesting feature of this calculation 
is that the second order correction to the eigenfunction is the term which describes the two gluon exchange in the 
t-channel such that the two gluons are in the octet. The fact that this term reggeizes precisely in the same way as 
the one gluon exchange is the essence of the celebrated bootstrap feature in high energy QCD 15[ . 



Bootstrap of the symmetric adjoint 

The two t-channel gluons in the previous calculation are in the antisymmetric adjoint representation. It is easy to 
show that two gluons in the symmetric adjoint also reggeize (see second paper in [8|). To see this let us consider a 
similar calculation but take the matrix R to be in the fundamental representation. As before we take 

G% = [ ah u tr[T a R F (u)] (60) 
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where r a are generators of SU(N C ) in the fundamental representation. The action of Hklwmij on this state is 
claculated just like before using 



[Ji(x),R F (y)}=T a R F (y)S 2 (x-y); [J%(x), R F (y)} = R F (y)T a S 2 (x - y) 
Using of completeness relation of fundamental SU(N) generators 



the properties of fundamental generators 



tr(T a T b ) = ^S ab , tr(T a T b T c ) = ^(d abc + if abc ) 
and the representation of an adjoint unitary matrix in terms of fundamental matrices 

R${z) = 2tr \r a R F (z)T b R F {z) 



one can write 



H 



KLWMIjGp 



2vr 2 



The first term on the RHS starts in the order (S/Sp) 3 . Thus the eigenvalue equation to second order reads 
aN. 



2tt 2 



K 



uuz y^u 



(61) 



(62) 



(63) 



(64) 



z <\> u \tr[\ - R F {z)R F {u)]tr[T a R F {z)] + N c tr [r a (R F (u) - R F (z))]j (65) 



J_ M , 1 ( if abc , d a6c| f ( S ^ * S 

5pl Sp a J 2 XJ ' J x - >y -\5p b (u,x-)6p c (u,y-) Sp b (z,x-) 6p c (z,y-) 



\^ + ^ ifabc+dabc U 

Ju \_ u Hu * Jx~>y- 



Sp b (u,x~) 5p c (u,y ) 



(66) 



Again the plane wave eq.flU) is the solution of this equation with the eigenvalue eq.(05]). The second order term 
proportional to the d abc tensor corresponds to exchange of two t - channel gluons in the symmetric adjoint represen- 
tation. As we have seen in the previous subsection, the antisymmetric octet (the f abc term) reggeizes. Thus eq. (|66[) 
tells us that the symmetric adjoint reggeizes by itself. This is another example of bootstrap at work. 



The bootstrap in the KLWMIJ/ JIMWLK approach. 



Since the bootstrap plays such an important role in the discussions of high energy amplitudes, it is worth while 
explaining how and why the bootstrap condition in the KLWMIJ approach is satisfied by fiat. 

First off, we note that the bootstrap condition, which leads to rcggeization of the two gluon exchange can be stated 
as the relation between the real part of the kernel of the BFKL equation and the reggeized gluon traiectorv|16j. 



w(fci) - w(k 2 ) - + fc 2 ) 



1 



K(ki,k 2 ,qi,q 2 ) 



(67) 



Here K is the real part of the BFKL kernel which evolves the color singlet two i-channel gluon state. It arises from 
the first term in Hbfkl eq.([28p. Expressing it in coordinate space in terms of the kernel K (x, y, z) we have 



K{xy, uv) = 2 



K uvz 5(x - u)8(y - v) - K uvy 5(x - u) - K uvx S(y - v) + K uvx S(x - y) 



(68) 



According to eq. (|43[) and eq. (|44p the Fourier transform of ui(k) into coordinate space is f3(x). Fourier transforming 
the bootstrap condition into the coordinate space we have 



P(x)6{y) - P(y)6(x) - (i(x)8{x - y) 



K ozS(x)6(y) - K 00v 5(x) - K 00x 5(y) + K 00x 5(x - y) 



(69) 
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We can derive this condition directly by considering the eigenvalue equation in the symmetric adjoint channel. Take 
the trial function in the form 



J uv u ru u rv 



(70) 



Acting on it by Hbfkl we derive the eigenvalue equation 

_ NcfX 
4tt 2 



2K uvz V(u, v) - 2K uzv ^(u, z) - 2K zvu V(z, v) + 2K zx J(u - v)V(z, x)\ + [/9 ao *(z, u) + P zu $(z, v) 

= u q *(u,v) (71) 

Assuming that the solution identical to the single gluon exchange 

*,(«,») = 5 2 {u-vy qu (72) 

exists for all and taking the integral over q we arrive at the configuration space condition eq. (|69| [22I ] . 

With f3(x) defined in eq. (f29| . this condition is clearly satisfied. Clearly, the bootstrap condition would be satisfied 
in the KLWMIJ approach for any functional form of kernel K xyz , since the relation between /3 and K, eq. (|29p is 
immutable. It appears simply due to normal ordering of Hbfkl written in the original form eq.(|26p. Thus even if 
K is modified in cq.(f26)). the bootstrap condition will still be automatically satisfied. One can contemplate several 



reasons for such a modification. First, higher order corrections in a s certainly lead to modification of K 17]. Another 
reason to consider a modification of K is the unphysical infrared behavior of perturbative gluon emission which leads 
to violation of Froissart bound 18]. Cutting off long distance tails of the Weiszacker- Williams field in the emission 
kernel if is a possible " phenomenological" solution of this problem 19]. One could question in principle the starting 
point of our discussion - eg. (1261) . However the form eg. (|26p is dictated by the hermiticity of Hbfkl- The real emission 
part is given by the first term in eq. (|28j) . By itself this term is not hermitian, and only with the gluon trajectory 
term, the second term in eq. ([28)) . the hermiticity of the Hamiltonian is restored. 

Thus we conclude that the bootstrap condition within the KLWMIJ framework is tantamount to the condition of 
hermiticity of the Hamiltonian which generates the rapidity evolution of the scattering amplitude. 

The Hermiticity of the Hamiltonian is related to the unitarity of the high energy evolution, even though the evolution 
equation is not a Schroedinger equation, but rather a diffusion type equation. The evolution Hamiltonian acts on the 
probability distribution W eq.Q. Since W[p] has the meaning of probability density, it must be positive definite. The 
eigenvalues of the evolution thus better be real, otherwise the "probability density" will develop an imaginary part 
even if one starts with a real and positive distribution at initial rapidity. This is assured if the evolution Hamiltonian 
is Hermitian. 

The origin of reggeization, including the gluon rerggeization, is in the i-channel unitarity. On the other hand the 
JIMWLK/KLWMIJ approach is formulated in s-channel and has the most natural interpretation as the evolution of 
the s-channel wave function. It is therefore interesting to see how t and s channel pictures are interrelated on the 
example of the gluon reggeization. The two pictures lead to the equivalent description, if the evolution in rapidity 
can be described by a hermitian Hamiltonian. 

We note that the bootstrap equation was used in [16[ to find the generalization of the gluon reggeization for the 
running QCD coupling case. The JIMWLK/KLWMIJ Hamiltonian is of course modified when the running is taken 
into account. The practical conclusion from our discussion in this subsection is the following: if we know how to 
include the running of the QCD coupling in /3 of Eq. (|29|1 . then Eq. (|28| can be used to generalize the full BFKL 
Hamiltonian to the running coupling case. This idea has been explored in [l6| and led to the "triumvirate" structure 
[201 ] . The same form of the Hamiltonian was derived recently in [20j by direct summation of the Feyman diagrams in 
the dipole approximation. Eq. (|28|) with 



= J_ [* a((( 9 -fc) 2 )a^) q* 

2ttA a(qi) &(q-k)* 



(73) 



gives the correct generalization of the JIMWLK/KLWMIJ Hamiltotian for running QCD coupling both for linear and 
non-linear term at any value of N c . 



SCREENING CORRECTIONS 



As we have seen above, expansion in powers of S/Sp is the expansion in number of gluons exchanged in the i-channel. 
For example the linear approximation of ea. (|42[) allows only one gluon exchange between the partons of the projectile 
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FIG. 2: KLWMIJ evolution in the approximation where only one gluon exchange in the f-channel is allowed. 

and the target. Diagrammatically the approximate diagonalization discussed in the previous section corresponds to 
summing the diagrams of Figj2] The evolution therefore allows for emission of an arbitrary number of gluons in the 
wave function of the projectile, but only for a single gluon exchange between the evolved projectile and the target. 
The terms quadratic in 5/Sp are represented in Fig|3l The vanishing of the last diagram on Fig|3]leaves the remaining 
contributions local in transverse coordinate and thereby ensures the reggeization of the two gluon exchange. The third 
order terms are not local anymore. In particular one encounters the diagrams of FigU which give a non vanishing 
bilocal contribution. 

Since the quadratic terms in the expansion of the wave function reggeize in the same way as the linear term, the 
significant corrections in a sense start form the cubic order. There is no reggeization of the third order terms in eq.(|54p 
nor eq. (|65|) . which is another way of stating that three gluon exchange is sensitive to screening corections. It is thus 
interesting to calculate the terms of order (S/Sp) 3 in the wave function. Although this can be done, the calculation is 
somewhat tedious. In this section we will perform a similar calculation, but the one that is easier implemented and 
has a somewhat more direct meaning in the framework of Hklwmij- Instead of expanding the wave function ^[R] 
in powers of S/Sp we will expand it in powers of R — 1. This is similar in spirit to [6(. Since the complete KLWMIJ 
wave function must depend on i?, this type of expansion is more direct. To leading order the expansion in R — 1 and 
S/Sp are equivalent. Beyond the leading order however they differ both on the calculational and conceptual levels. 

Expansion in powers of R— 1 corresponds to the expansion in the number of the projectile partons which participate 
in scattering. To linear order in R — 1 our approximation allows only one parton in the projectile wave function to 
scatter off the target, but it can scatter by exchanging an arbitrary number of i-channel gluons. The diagrams which 
are resummed in this approximation are depicted on Fig. 4. Clearly the leading order expansion in S/Sp is subsumed 
in the leading order expansion in R — 1, since if we allow only one i-channel gluon, we also allow only one parton of 
the projectile to scatter. In higher orders it is not the case anymore. In this section we will calculate 0[(R — l) 2 ] 
correction to the wave function of eq. (|60|) . 

Recall eq.(|65l> 

H KLW MU J <p u tr[T a R F (u)} = £^ J K uuz <j) u {«r[l - R F {z)R F {u)]tr[T a R F (z)} + N c tr [r a (R F (u) - R F (z))]} 
= 2^2 / KuuzK i Nctr [t ° (i?F(u) ~ + M^H") - l]tr[T a R F (z)} - tr[R F {z) - l]tr[r a R F (z)} + ...}(74) 

The linear term on the RHS is the familiar reggeization term. We see thus that reggeization is the property of an 
arbitrary number of gluon exchanges, as long as all the ^-channel gluons couple to the same parton (in this case 
quark) in the projectile. The second term on RHS is second order in R — 1. Note that in terms of the expansion in 
S/Sp it actually starts with the cubic order. Clearly, to correct the eigenfunction we have to add to our original G F a 
term of second order in R — 1. Guided by the form of the RHS of eq. (|T4")) we take the wave function to second order 




FIG. 4: A nonlocal in transverse plain contribution to the eigenfunction with three t channel gluon exchanges, 
of the form 

Gp = G\ + Gl (75) 
G\ = \ 4, u tr\T a R F {u)\ 

J U 

G a 2 = f ^(u,v)tr[R F (v)-l]tr[T a R F (u)} + ^- f i>(u, v)tr[R F (v) - l\tr[ T a R F {u)] 

We will see that this ansatz is general enough to satisfy the eigenvalue equation to second order in the large N c limit. 
In the rest of this section all the matrixes R are in fundamental representation and we drop the subscript F for 
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FIG. 5: Reggeization diagramms with arbitrary number of £-charmel gluons that couple to the same parton in the projectile. 



convenience. 

The action of Hklwmij on G is straightforward to calculate. After some algebra we find 

Hklwmij j i/>(u, v)tr[R\ - l]tr[ T a R u ] = 
= i>(u,v)I^K uvz {tr[RlRlR z T a R u ]+tr[RlR u T a ^ 

-K vvz [tr{R\)tr[{Rl - l)(R z - 1)] + tr{R\ - l)tr{Rl - 1) 
+tr(R\ - l)tr(R\ - 1) + N c [tr{R\ - 1) + tr(R z - l)]}tr(r a i? u ) 
-K uuz \tr[R\ - l]tr[{R\ - 1){R U - 1)] + tr[R\ - l)tr[R\ - 1] 
+tr[R\ - l]tr[R u - 1] + N c tr[Rl - l]}tr[r a i? M ] 

+N c K uuz tr[Rl - l]tr[T a R u \ 1 (76) 

Keeping only [0(R — l) 2 ] terms and taking large N c limit for simplicity, we have 

Hklwmij J ^(u,v)tr[Rl-l}tr[T a R u }= (77) 

= 2^/ N c\^{u,z)K vvz tr[l- R v ]tr[r a R u ]+ \^{u, z)K vvz + ^{z, v)K uuz - ^(u, v)K uuz ] tr[l - Rt]tr[ T a R u ] j 
Similarly to quadratic order in the large N c limit 

Hklwmij J tr[R v - l]tr[T a R u ] = (78) 
= ^2 / Nc fa u ' z ) K wztr[l - Rt]tr[T a R u ] + (${u, z)K vvz + $(z, v)K uuz - $(u, v)K uuz ^j tr[l - R v ]tr[T a R u ] } 
Finally we have 

HG a 2 = 7^1 {^{u 7 z)K vvz +^{u,z)K vvz +i:{z,v)K uuz -^(u,v)K uuz ^tr[l- Rl]tr[T a R u ] 



+ 



2tt 2 



/ j -0(u, z)K vvz + ip(u, z)K vvz + %b(z, v)K uuz - 4>(u, v)K uuz \ tr[l - R v }tr{r a R u } (79) 



We determine the functions ip and ip by solving the eigenvalue equation 

HklwmijG f = loqG f 



(80) 
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with the eigenvalue ujq given by eq. (|45l) . As we have seen in ea (|74p the linear terms reggeize and cancel between the 
left and right hand side. For the quadratic terms we are then left with the equation 







a 






2^2 
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Q 
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2^2 






W Q 






N c 


4>u = 


exp{iQu} 


a 
2^2 


I 


ip(u, 


a 

2^ 


I 


4>(u, 



K uuz <j> u [tr{l - R\) + tr(l - R u )]tr(r a R z ) 

Z 

| ip(% z)K vvz + $(u, z)K vvz + Tp(z, v)K uuz - ip(u, v)K uuz > tr[l - Rl}tr[r a R u ] 

v,z 

I ip(u, z)K vvz + ip(u, z)K vvz + ip(z, v)K uuz - $(u, v)K uuz \ tr[l - R v }tr[r a R u ] 

v.z 

[i/j(u,v)tr(Rt - l)tr{T a R u ) + ^(u,v)tr(R v - l)tr(T a R u )} 



(81) 



a 
2^2 

a 
2^2 J 



ip(u, v)K uuz - %j){u, z)K vvz - tj}(z, v)K uuz - ip(u, z)K vvz - -^-^>{u, v) 
ip(u, v)K uuz - i/)(u, z)K vvz - t/j(z, v)K uuz - ip(u, z)K vvz - -^-4>(u, v) = -^K vvu <f) v 



K zzu (f> z S(u - v) (82) 



These equations are easily solved in momentum space. Defining 

V>(u,u) = J d 2 kd 2 qe l( - qu+kv ^(q,k); ^{q,k) = J d 2 ud 2 ve-^ qu+kv ^{u,v) 
it is starightforward to Fourier transform eq. ([82} . For example 

s- 4(<zu+fct,) I K vvz ^(u, z) = In ^(q, k); etc. 



(83) 



(84) 



(85) 



The ultraviolet cutoff A is needed to regularize the divergence in the Fourier transfrom of l/x 2 . Similarly Fourier 
transforming the rest of the terms we get 



A 2 



A 2 , A 2 



Q 



k 2 

A 2 A 2 A 2 \ - A 2 A 2 



A 2 



A 2 



ln 7^2 - ln ~2 ~ ln tt ^te fc ) - ln k ) = ln ty? s (Q-q-k) 



k 2 
k 2 



Q 



(86) 
(87) 



The solution is found straightforwardly as 



In ^ -In 



In^+ln^ 

i/>{q,k) = — V2 — ^— ^ — . 2 8 (Q - q - k) -^a^oo S (Q-q-k) 

In qj - In -r - 2 ln w 

Transforming back to the configuration space we find the eigenfunction to second order in R — 1 and in the large N c 
limit 



G F = I e lQu tr[T a R F (u)) 1 - — tr[R F (u) - 1] 



1 



(89) 



CONCLUSIONS 



To summarize we have discussed relationship between the KLWMIJ Hamiltonian and its BFKL limit. Eigenfunc- 
tions of H klwmij when expanded to leading order in S/Sp become eigenfunctions of Rbfkl- It is however difficult 



15 



to determine which eigenfunctions of Hbfkl become normalizable eigenfunctions of Hklwmij when the Taylor 
series is resummed to all orders. 

The relation of course pertains only to functions (functionals) which are expandable in Taylor series in 5/Sp. We 
know from the general discussion of |lOj that Hklwmij also has eigenfunctions which do not have such an expansion. 
Those are states close to the black disk limit. For those states the pertinent expansion is in powers of p. We note 
that Hb fkl is self dual under the transformation 



It thus contains eigenstates whose eigenfunctions are monomials in p rather than 6/Sp. The self duality ensures that 
those have exactly the same spectrum as the eigenfunctions we discussed in the bulk of this paper, even though these 
states are indeed formally close to the black disk limit 

The duality transformation eq. (|5U]) is the linearized version of eq. (fT5|) which transforms Hklwmij to Hjimwlk ■ 
Thus the same relation as discussed above exists between the second set of the eigenfunctions of Hbfkl and the 
eigenfunctions of Hjimwlk- 

We have also discussed the gluon reggeization and the appearance of the bootstrap condition in the KWLMIJ 
formalism. The bootstrap condition is direct consequence of Hermiticity of Hkwlmij and as such is a necessary 
attribute of the approach. Any modification of the emission kernel K xyz does not ruin the bootstrap property. Since 
the JIMWLK picture is a s-channel one while the reggeization stems from the t channel unitarity, we conclude that 
the Hermiticity of Hkwlmij is the property that reconciles these two approaches. 

Further we have discussed expansion of the eigenfunctions in powers of R — 1 rather than powers of S/Sp. This 
corresponds to expansion in the number of partons in the projectile wave function which participate in the scattering. 
We have calculated 0[(R — l) 2 ] correction to the reggeized gluon wave function and found that is has a very simple 
form in the large N c limit. 

We note in this respect that while the eigenvalue determines the rapidity dependence of the scattering amplitude, 
the functional form of the wave function determines the impact factor. This follows from the expansion eq.© which 
can be written as the overlap of the "wave function" characterizing the projectile hadron at the initial rapidity and 
the eigenfunction of the evolution Hamiltonian 



Thus for example the eigenfunction eq. (l5"4l will not have zero overlap with a single quark state, while that of eq. (j6"0")) 
will not overlap with a gluon projectile, since their incoming color representations cannot be combined into a singlet. 
Thus in order for the amplitudes of both, quark and gluon projectiles to have the same energy dependence, there must 
be a degeneracy in the spectrum of Hklwmij- We indeed saw this degeneracy explicitly since both eigenfunctions 
eq. (l54l and eq. (|60|) correspond to the same eigenvalue. On the formal level this degeneracy is the consequence of the 
fact that perturbative vacuum breaks the global symmetry of Hklwmij as SUl(N c ) ® SUr(N c ) — > SUv(N c ), the 
reggeized gluon being " the Goldstone boson" associated with this breaking Q . 
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In order to derive the BFKL Hamiltonian we have to expand R ab {x) and Jl(r)( x ) m powers of Spa( l x - y The 
expansion of R is straightforward. To expand Jr and Jl we will use their commutation relations with R(x). We start 
with the KLWMIJ Hamiltonian 




(90) 



7i = <WV |¥.) 



(91) 
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Expansion of R is straightforward. To second order we have 



Rf = 5»P + Tf / du- — + T^Tf / du x / du 2 

Jo dp c (u,u 1 ) Jo Jo 



(93) 



2 Sp c (u,u 1 ) 6p d (u,u 2 ) 
To expand Jl we use the commutation relation 

[Jl(x),R(y)} = T a R(y)S 2 (x-y) (94) 
It is easy to check that to second order in 5/Sp the following expression satisfies the correct commutation relation 

5 



JIM 



du l p a (u, u l ) - T* K \ du x I du 2 p x {u, u 2 ) 



o Jo 

b 7b , 







5p K (u, u 1 



Now using the fact that J£(w) = R^f Jr{u) we have also 



J a R (u) = - du x p a (u,u x ) + T^ du x du 2 pX(u,u 2 ) 



5p K (u, u 1 ) 



(95) 



(96) 



Now the expansion of the amplitude Q reads 



(x 


- z)i 
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-zf 
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- Z)i 


[x 


-zf 
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- z)i 


(x 


-zf 



[R ab (z)-R ab (x)}j b R (x) 

rpOtfl 



1 f 1 6 

dx~ dx 2 p a (x,x-)—— , „., , , 

i) Jo opP{x,x 2 ) Jo Jo 

5 S 



1 - r 1 

dxi I dz x p a {x, x x ) 



Opx Opt 



(97) 



with f Q dx p a (x,x ) = and J Q dx Sf>a ^ x = jj^; ■ Finally we can write the BFKL Hamiltonian as in ea(|26p. 
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